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Abstract 

We prove that the weakly singular, non-linear convolution integral 
equation \x — y\~^f{y)dy = where 0 < A < n, and p = 

2nl{2n — A) has at least two non-equivalent solutions. This answers a 
problem of Elliott Lieb. We also prove certain orthogonality relations 
among linear differential forms with constant coefficients related to the 
corresponding type of convolution operators. Finally, we discuss the 
regularity of the solutions of such non-linear integral equations over 
not necessarily bounded open subsets of R". 


1 Introduction of the results 

We divide the results in this paper into three parts. In the first part we 
note that a certain non-linear integral equation introduced by Elliott Lieb 
has at least two essentially different solutions. One of the solutions has an 
isolated singular point where the function tends to infinity. 

The second part presents a multitude of integral identities each con¬ 
necting two solutions of Lieb’s equation. These integral identities contain 
the images of the two solutions under linear differential forms with constant 
coefficients. They point out to two properties connecting any two such so¬ 
lutions. One property is an orthogonality property and the second property 
is a certain commutativity. 

Finally, in the third part we present regularity results of the solutions 
to the Lieb integral equation. The bottom line is that except for isolated 
singularities were the solutions tend to infinity, they are, in fact smooth 
functions elsewhere. One can compute effectively their smoothness degree. 

(I) In d] Elliott H. Lieb computes the sharp constants for certain param¬ 
eter values in the Hardy-Littlewood-Sobolev inequality. The paper also deals 
with related inequalities: the Sobolev inequality, doubly weighted Hardy- 
Littlewood-Sobolev inequality and the weighted Young inequality (as A. 
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Sokal called it). Theorem 3.1 on page 359 computes the unique maximiz¬ 
ing function and the sharp constant for certain values of the parameters of 
the first inequality. The maximizing functions should satisfy the following 
integral equation 

/ \x-y\~^ fiy)dy = f{x)P-^, 0<A<n, p= ‘^^ ( 1 . 1 ) 

2 n — a 

However, as claimed on page 361 of [T] ”We do not that (3.9)” (the above 
(jl.ip equation) ’’has an (essentially) unique solution-even if we restrict to 
the SSD category-and we shall offer no proof of this kind of uniqueness. 
This is an open problemP . Indeed there is no uniqueness for we have the 
following: 

Theorem 1.1. The following function is a solution of equation (EH).- 

/(x) = C(n,A)|x|-(-V2), (1.2) 


where: 

C{n,\)= (1.3) 



Remark 1.2. Lieb proved that his equation O) has the following solution 
(which is unique among the maximizing functions of the inequality he was 
considering): 


hix) = L{n, A)(l + |x|2)—/P = L(n, A)(l + |x|2)-(-V2). ( 1 . 4 ) 

Here L{n, A) is a constant depending on the dimension and on the parameter 
A. We note that the solution in Theorem 1.1 is singular at the origin and is 
not a reflection of Lieb’s solution. It is certainly not a conformal image of 
it {fr is non-singular and bounded). Thus Lieb’s equation (II.ip exhibits at 
least two non-equivalent solutions. 

Based on the coming results we can deduce many integral identities. 
Here is an example: 

Corollary 1.3. 

L{n, f |^|-(n-A/ 2 )^^ ^ \xf)-^/^dx = 

(1.5) 

= L(n,A)C(n,A)^/(2n-A) f |a;|-V2(i ^ |^|2)-(n-A/2)^^^ 


2 



We note the symmetric relations between the two solutions within the 
last identity. It is a kind of commutativity that takes the power p — 1 into 
an account. 

(11) The identity in the last corollary is the zero’th integral identity out 
of infinitely many possible integral identities that connect two solutions of 
the Lieb integral equation (ini). 

Definition 1.4. For a = (oi,..., an) G (2^ U {0})”' we denote: 



where Oq, G M. 


Here are some integral relations between solutions of equation 0. 
Obvious generalizations hold true for other similar kernels. These include 


some orthogonality relations and some commutativity relations: 

Theorem 1.5. Let f{x) and g{x) he solutions of equation Then 

assuming the convergence of the integrals below, \/a,l3 G (Z+ U {0})"' we 
have: 




(1.7) 



[ Df\f{x))D(-\f{x)P-^)dx= [ D^^\f{x))D<f\f{x)P-^)d: 

JR^ JR" 


\x = 0, 


( 1 . 8 ) 


So it is natural to make the following: 


Definition 1.6. 



O = {^ b0Df\-) \\//3eJ, (-1)1^1 = -1, b0 G M}. 


/3gJ 
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Theorem 1.7. Let f{x) and g{x) he solutions of equation 111.1]) . Suppose 
that A = Ae + Ao, 12 = Qe + where Ag, lig € E, and Aq, 12o G O. Then: 

(a) Assuming convergence of the integrals below we have: 

(1.9) 

/ Kf)^{f~^)dx = [ A{fP~^)Q{g)dx, 

JR^ JR" 

[ Ae{f{x))Ao{f{x)P~^)dx = [ Ae{f{xy~^)Ao{f{x))dx = 0. 

JR" JR" 

(b) Assuming the convergence of the integrals below we have: 

(1.10) 

[ A{f)n{fP-^)dx = [ A{fP-^)Q{f)dx = 

JR" JR" 

[ Aeif)ileiF~^)dx + f Aoif)noif^~^)dx = 

JR" JR" 

= / AeiF-^)n,if)dx + [ Ao{fp-^)noif)dx. 

JR" JR" 

(III) We turn our attention to the smoothness of the solutions of Lieb’s 
equation dni). We recall some notations and definitions from [3]. Let 
G C be an open and bounded set. For a A G G. Vainniko introduces 
a weight function: 

( 1 for A < 0 

rcA(j;) = s (1 + I logp(x)|)“^ for A = 0 , x € G, 

[ p{x)^ for A > 0 

where p{x) = mfy^gc \x — y\ is the distance from x to the boundary dG of 
G. Let m G U {0}, i/ G E satisfy u < n. We define the space C'™’’^(G) as 
the collection of all m times continuously differentiable functions u : G —)■ E 
(or C) such that: 

\\u\\m,u = ^ SUp{w\c,\-(n-u){x)\DaU{x)\) < OO. 

I 1^ x£G 
\a\<m 

So C^''^{G) contains all the m times continuously differentiable functions u 
on G whose derivatives near the boundary dG can be estimated as follows: 

{ 1 for |q;| < n — V 

1 + I log/9(a:) I for |q;| = n — , x G G, \a\ < m. 

p(a:)”“^“l“l for |q;| > n — 
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Remark 1.8. 1) The function || • \ \m,i> on the space C'^’^{G) is a norm. 

2) The space || • \\m,u) is complete, i.e. it is a Banach space. 

Consider the following integral equation: 

u{x) = / K{x,y,u{y))dy + f{x), xeG. (1.11) 

Jg 

we assume that the kernel K{x,y,u) is m times (m > 1) continuously dif¬ 
ferentiable with respect to x,y,u iov x € G, y € G, x y, u G M. We 
also assume that there is a real number v < n, such that, V fc € and 
a, P G (Z'*')”' for which A: -|- |a| -|- |/3| < m, the following inequalities hold: 

( 1 . 12 ) 

f)k f 1 for 1 / + |a| < 0 

D^D^^'^^-^K{x,y,u) <bi{u) < 1 + | log |x - ?/|| for |a| = 0 , 

^ I |x —for + |a| > 0 

(1.13) 

f)k ok 

{ 1 for z^ + |a| < 0 

1 I log |x - y|| for v + \a\=Q . 

|x —for v + \a\>t) 

The functions 6 i : M —)• and 62 : are assumed to be bounded on 

every bounded region of M and respectively. The notation of G. Vainikko 
reads as follows: 

(1.14) 

We can now conveniently quote the result we need from [3]: 

Theorem 8.1.((3j) Let f € G^’^{G) (in lll.lll) } and let the kernel K{x, y, u) 
satisfy inequalities and 11.13)) . If the integral equation 11.111) has a 

solution u G L°°{G), then u € C'™'’^(G). 
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Remark 1.9. 1) There is a companion result (Theorem 8.2) in [3], but we 
will not use it here. 

2) We note that if we restrict Lieb’s integral equation (11.111 to a bounded 
open G C M", then it satisfies the assumptions of Theorem 8.1. in [3]. In 
this case f{x) = 0 € C°°{G) and K{x, y, u) = \x — y\~^ is independent of u. 
Also we note that: 



and so in inequalities (11.1211 and (|1.13ll the only interesting values of (/c, |q;|, |/3|) 
are k = 0, |/3| = 0, |a| < m and within the inequalities we start with 
u + \a\ = \. Thus we can make effective calculations of the smoothness de¬ 
gree of a bounded solution of the restriction to G of Leib’s integral equation 


m- 


3) We note that any solution f{x) of the Lieb integral equation (11.111 . must 
decay to zero at infinity, in order for the integral 



to converge at infinity. Thus any solution is bounded outside a ball Bn{R) = 
{x € M" I |x| < i?} for a large enough radius R. Thus such a solution can 
have singularities only within the ball Bn{R), and the solution tends to 
infinity at each such a singularity, otherwise by Theorem 8.1 in [3] it will be 
smooth (of some degree) at such a point. 

We thus have the following: 

Theorem 1.10. Let f{x) he a solution of Lieb’s integral equation M.l\) . 
Then lim| 3 ,|^Qo f{x) = 0, and there is a ball Bn{R) such that f{x) is bounded 
for X ^ Bn{R), f{x) can have a finite set of singularities inside Bn{R), and 
f{x) tends to infinity when x ^ s for each singular point s of f{x). 


□ 


2 Lieb’s integral equation has at least two non¬ 
equivalent symmetric decreasing solutions, The¬ 
orem 1.1 

A proof. We use the following formula for the Fourier transform, [2]: 


|y|-" 


Ijyr exp(-2«x ■y)dy= (| - 0 /r (0 } 


X 
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where 0 < u < n. We compute the Fourier transform of our integral: 


([ \t-y\~^f{y)dy) (x) = ( [ \t-y\-^{C{n,\)\y\ (x) = 

= C{n, X)( [ lyl"^ exp {ix ■ y) dy^ ( [ exp (ix ■ y) dy \ = 

Vjr" / \Jmx J 



= C(n,A)^/(2n-A) j'|y|-(n-V^/(2n-A)^ (3.) = 


= (’{C(n.A)|!,|-(”-V2))V(2«-A)^ ^ 

= m = (/fer') M- 

Hence: 

[ \x-y\~^f{y)dy = f{xf-^. 

JR’^ 

□ 


3 Integral relations between two solutions of equa¬ 
tion (ll.ll) . orthogonality and commutativity, Corol¬ 
lary 1.3, Theorem 1.5 and Theorem 1.7 

In this section we are interested in the form of integral relations between 
two solutions of the Lieb integral equations. Thus we will not bother with 
convergence issues and just formally expand the formulas in Corollary 1.3, 
in Theorem 1.5 and in Theorem 1.7. We note that the identity in Corollary 
1.3 follows by Theorem 1.1 and by the case a = /3 = 0 in equation (II.6p of 
Theorem 1.5. Also Theorem 1.7 follows by Theorem 1.5. Hence we need to 
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prove only Theorem 1.5: 

We start with Lieb’s integral equation, (|l.ll) and perform on it a partial dif¬ 
ferentiation with respect to Xj. We formally differentiate under the integral 
sign. The result is: 


/ ^ fiy)dy 

J]^n OXj V / 


df{xy ^ 

dxn 


We note that: 


Hence: 


dx 


^ (k - !/r") =- !/r"). 




Jmp 9y 

By integration by parts we deduce the following: 


dxi 


/ k - 2/1 


-xdfiy), df{x)P 
dy = 


-1 


dyj 


dxi 


We iterate this argument and obtain: 


dykdyj ^ dxkdxj 
Now an inductive argument implies: 

[ \x- y\~^ Dy'> (/(y)) dy = {f{xy~^) . 

Another outer induction gives, finally: 

[ \x-y\~^ A{f{y))dy = A{f{xy~^) . 

JR" 

Next, let g{x) be one more solution of equation (II.ip . i.e.: 

[ \x-y\~^ g{y)dy = g{xf-^. 

Jr" 

By what we have already done, we have: 

f \x- y\~^ {g{y)) dy = . 


(3.15) 


(3.16) 
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We now use the double integration technique. We multiply equation ()3.15l) 
by D^^\g{x)) and integrate ... dx: 

[ iaix)) [ \x-y\~^D^y^f{y))dydx= [ Df\g{x)) D^^^f{x)P-^) dx. 

Jr^ Jr^ 

Reversing the order of integration on the left hand side gives: 

[ D^aHfiy)) [ \x-y\~^Df\g{x))dxdy=[ D^y^f{y))D^^\g{y)P-'^)dy. 

JR^ JR^ Jr^ 

Changing on the right hand side the name of the variable from y to x, gives: 

[ Df{g{x))D^^^[f{x)P-^)dx= I D^^){f{x))Df [g{x)P-^)dx. 

This proves equation (|1.6p . This is the commutativity part. We now prove 
orthogonality. We start with equation (|1.6p in the special case f{x) = g{x). 

We get: 

[ if{x)) {fix)P-^) dx= I (/(x)) (/(xf dx. 

jR" JR" 

Integration by parts gives: 

[ Dl^Hfix))D^^'> {f{x)P-^)dx = {-l)\‘^\ f fix) ■ D^^lp{f{x))dx, 

JR" JR" 

while 

[ {fix)P~^) (fix)) dx = (- 1 ) 1^1 f fix) • T>2/3 

JR" JR" 

We proved equation (11.71) . Equation (|1.8I) is a consequence of equation (11.61) 
and equation ()1.7I) . This proves Corollary 1.3, Theorem 1.5 and Theorem 
1.7. □ 
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